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Using a high resolution spectrum of the secondary star in the neutron star binary Cen X-4, 
we have derived the stellar parameters and veiling caused by the accretion disk in a consistent 
way. We have used a \ 2 minimization procedure to explore a grid of 1 500 000 LTE synthetic 
spectra computed for a plausible range of both stellar and veiling parameters. Adopting the best 
model parameters found, we have determined atmospheric abundances of Fe, Ca, Ti, Ni and Al. 
These element abundances are super solar ([Fe/H] = 0.23 ± 0.10), but only the abundance of 
Ti and Ni appear to be moderately enhanced (> la) as compared with the average values of 
stars of similar iron content. These element abundances can be explained if the secondary star 
captured a significant amount of matter ejected from a spherically symmetric supernova explosion 
of a 4 Mq He core progenitor and assuming solar abundances as primordial abundances in the 
secondary star. The kinematic properties of the system indicate that the neutron star received 
a natal kick velocity through an aspherical SN and/or an asymmetric neutrino emission. The 
former scenario might be ruled out since our model computations cannot produce acceptable fits 
to the observed abundances. We have also examined whether this system could have formed in 
the Galactic halo, and our simulations show that this possibility seems unlikely. We also report 
a new determination of the Li abundance consistent with previous studies that is unusually high 
and close to the cosmic Li abundance in the Galactic disk. 
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ABSTRACT 



l. 



Introduction 



"Based on observations obtained with UVES at VLT 
Kueyen 8.2 m telescope in programme 65.H-0447 

1 Consejo Superior de Investigaciones Cienti'ficas, SPAIN 



The optical counterpart of the low mass X- 
ray binary (LMXB) Cen X-4 (V822 Cen) was 
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Fig. 1. — Observed spectrum of the secondary star of Cen X-4 (top) and of a properly broadened template 
(HD 216803, bottom). 



firstly identified during an X-ray outburst in f 979 
(Kaluzicnski et al. 1980), ten years after the X-ray 
outburst detected by the satellite Vela 5B during 
an X-ray outburst in 1969 (Conner et al. 1969). At 
that moment, it resembled an my ~ 13 magnitude 
blue star that had increased in brightness by ^6 
magnitudes from its pre-outburst value (Canizares 
et al. 1980). About three weeks after the 1979 out- 
burst a type I X-ray burst was observed, which 
indicated that the X-ray source is a neutron star 
(Matsuoka et al. 1980). Two weeks later the op- 
tical counterpart had faded to my = 18.2 mag 
and spectral observations revealed the presence of 
a late-type companion star (K3-K7) in the sys- 
tem (van Paradijs et al. 1980). Later optical stud- 
ies determined the orbital period, P = 0.629 d 
and a system center-of-mass velocity with respect 
to the Sun of 137 ± 17 km s" 1 (McClintock & 
Remillard 1990). The corresponding mass func- 
tion was then estimated at f(M) = 0.23 M , 



which further supported the argument for a neu- 
tron star primary (rather than a black hole, Mc- 
Clintock & Remillard 1990). More recently, from 
measurements of the orbital inclination and binary 
mass ratio, the compact object mass has been es- 
timated at M NS = 0.5-2.1 M Q (Shahbaz et al. 
1993) and the companion star mass, in the range 
0.04 < M 2 < 0.58 M Q (Torres et al. 2002). 

The evolutionary status of this system has 
been widely discussed in the literature. Whereas 
Chevalier et al. (1989) and McClintock & Remil- 
lard (1990) had suggested a scenario in which the 
secondary star was a "stripped giant" with a very 
low mass (M 2 ~ 0.1 M Q ), Shahbaz et al. (1993) 
favored the case for a subgiant. In both cases 
the secondary star would be an evolved star that 
had spent more than 80 per cent of its life (i.e., 
~ 10 9 -10 10 yr depending on its initial mass). On 
the contrary, the compact object in this system is 
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a neutron star, probably the remnant of the SN 
explosion of the initial massive component of the 
binary. The ejected material in the SN explosion 
is composed of heavy nuclei as products of the 
hydrostatic and explosive nucleosynthesis in the 
massive star. 

Recent studies of chemical abundances of sec- 
ondary stars in low mass X-ray binaries have 
opened up a new window of information on the 
later stages in the evolution of massive stars. A 
significant amount of the ejected matter in any 
supernova explosion that formed the compact ob- 
jects in these systems could be captured by the 
companion stars. These stars could therefore show 
anomalous abundances as signatures of the ex- 
plosion. Detailed chemical analysis of the at- 
mospheres of the secondary stars could constrain 
many parameters involved in supernova explosion 
models, such as the mass cut, the amount of fall- 
back matter, any possible mixing, and explosion 
energies and geometries. Indeed, the chemical 
analysis of the secondary stars provides evidence 
for a supernova event in the history of the sys- 
tems Nova Sco 94 (Israelian et al. 1999; Brown et 
al. 2000; Podsiadlowski et al. 2002) and A0620-00 
(Gonzalez Hernandez et al. 2004). 

In this paper we present the chemical abun- 
dance analysis of the secondary star of the LMXB 
Cen X-4 and discuss these results in the context 
of possible evolutionary scenarios. 

2. Observations 

We obtained 20 spectra of Cen X-4 with the 
UV- Visual Echelle Spectrograph (UVES) at the 
European Southern Observatory (ESO), Observa- 
torio Cerro Paranal, using the 8.2 m Very Large 
Telescope (VLT) on 2000 April 25 and 2000 June 
9, covering the spectral regions AA4800-5800 A 
and AA5800-6800 A at resolving power A/<5A ~ 
43, 000. Short exposures (718 s) were chosen in or- 
der to minimize possible smearing of spectral lines 
associated with the radial velocity change during 
its orbital motion, which we estimated to be less 
than 11 km s _1 . 

The spectra were reduced in a standard manner 
using the UVES reduction package within the MI- 
DAS environment. The radial velocity for each 
spectrum was obtained from the cphemcris re- 
ported in Casares et al. (2005, in preparation). 



The individual spectra were corrected for radial 
velocity and combined in order to improve the 
signal-to-noise ratio. After binning in wavelength 
in steps of 0.21 mA the final spectrum had a signal- 
to-noise ratio of 110 in the continuum. This spec- 
trum is shown in Figure 1. 

3. Chemical Analysis 
3.1. Stellar Parameters 

The contribution of the emission from the ac- 
cretion disk to the total flux in Cen X-4 has been 
estimated at ~ 25-30% and 40-50% in the V and 
B bands, respectively, whereas it is at least 80% 
at A3800 A(Chevalier et al. 1989; Torres et al. 
2002). As in Gonzalez Hernandez et al. (2004), we 
have obtained the veiling, together with the stel- 
lar atmospheric parameters, using synthetic spec- 
tral fits to the high resolution spectrum of the sec- 
ondary star in Cen X-4. First, moderately strong 
and relatively unblended lines of several elements 
of interest were identified in the high resolution 
solar flux atlas of Kurucz et al. (1984). We se- 
lected several spectral features containing in total 
29 absorption lines of Fe I with excitation poten- 
tials between 0.5 and 4.5 eV. In order to compute 
synthetic spectra for these features, we adopted 
the atomic line data from the Vienna Atomic Line 
Database (VALD, Piskunov et al. 1995) and used 
a grid of local thermodynamic equilibrium (LTE) 
models of atmospheres provided by Kurucz (1992, 
private communication). Synthetic spectra were 
then computed using the LTE code MOOG (Sne- 
den 1973). To minimize the effects associated with 
the errors in the transition probabilities of atomic 
lines, we adjusted the oscillator strengths, log gf 
values of the selected lines until we succeeded in 
reproducing the solar atlas of Kurucz et al. (1984) 
with solar abundances (Anders & Grevesse 1989). 

We generated a grid of synthetic spectra in 
terms of five free parameters, three to character- 
ize the star atmospheric model (effective tempe- 
rature, T e ff, surface gravity, logy, and metallic- 
ity, [Fe/H]), and two further parameters to take 
into account the effect of the accretion disk emis- 
sion on the stellar spectrum. This was assumed to 
be a linear function of wavelength and thus char- 
acterized by two parameters: veiling at 4500 A, 
/4500 = F^ 8 l°/F^ 5 c oa , and the slope, m . Note that 
the total flux is defined as F tota \ = F d i S k + 
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Table 1 

Ranges and steps of model parameters 



Parameter 


Range 


Step 


T e f{ 


4000 -> 5500 K 


100 K 


logg 


2.5 -> 5 


0.1 


[Fc/H] 


-1. -> 1 


0.05 


/4500 


0^2 


0.10 


ma 


-> -0.00082 


-0.00002 



Table 2 

Uncertainties in the abundances of the secondary star in Cen X-4 



Element [E/H] LTE A,, A T<)ff Ai og g A ve ii ing Atotal 



Al 


0.30 


0.12 


0.05 


0.05 


0.10 


0.17 


Ca 


0.21 


0.07 


0.10 


-0.12 


0.02 


0.17 


Ti 


0.40 


0.06 


0.13 


0.04 


0.09 


0.17 


Fc 


0.23 


0.09 


0.02 


0.02 


0.03 


0.10 


Ni 


0.35 


0.06 


0.06 


0.02 


0.05 


0.10 


Li* 


3.06 


0.12 


0.10 


0.02 


0.25 


0.29 



*Li abundance is expressed as: 

loge(Li) LT E - log [AT (Li) /AT(H)] lte + 12 



The total error was calculated using the following formula: 
Atotal = J K + A|T + Af og ~ + A*~ 



NOTE. — The errors from the dispersion of the best fits to different features, 
A CT , are estimated using the following formula: A CT — a/y^N, where a is the 
standard deviation of the measurements. Total errors also take into account 
the uncertainties associated with the stellar parameters and the veiling. 
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Fig. 2. — Distributions obtained for each parameter using Monte Carlo simulations. The labels at the top of 
each bin indicate the number of simulations consistent with the bin value. The total number of simulations 
was 1000. 



where i*disk and F scc are the flux contributions of 
the disk and the continuum of the secondary star, 
respectively. These five parameters were changed 
according to the steps and ranges given in Table 
1. A rotational broadening of 44 kms -1 , a limb- 
darkening e = 0.65 and a fixed value for the mi- 
croturbulence, £ = 2 kms -1 , was adopted. 

The observed spectrum was compared with 
each of the 1 500 000 synthetic spectra in the grid 
via a x 2 minimization procedure that provided the 
best model fit. Using a bootstrap Monte-Carlo 
method we defined the la confidence regions for 
the five free parameters and established the most 
likely values: T off = 4500±100 K, logs = 3.9±0.3, 
[Fe/H]= 0.4 ± 0.15, / 45 oo = 1-85 ± 0.10, and 
too = -0.00071 ± 0.00003. Confidence regions de- 
termined using 1000 realizations arc shown in Fig- 
ure 2. This spectroscopic determination of the sur- 
face gravity is consistent with derived surface grav- 
ities from the mass and radius of the secondary 
star that lie between logg ~ 3.6-4.0. 



3.2. Stellar Abundances 

Using the derived stellar parameters we analy- 
zed several spectral regions where we have iden- 
tified various lines of Fe, Ca, Al, Ti, Ni, and Li, 
often blended with Fe lines. Each of these spec- 
tral regions was carefully normalized using a late 
type star template (HD 216803) properly broad- 
ened with the rotational profile used for the sec- 
ondary star in Cen X-4 for comparison. We deter- 
mined the abundances of these elements by com- 
paring the observed spectrum with a grid of syn- 
thetic spectra through a \ 2 minimization proce- 
dure. For these spectral syntheses we modified 
the element abundances while the stellar parame- 
ters and the suitable veiling factor for each spec- 
tral region were kept fixed (for further details, see 
Gonzalez Hernandez et al. 2004). 

In Figs 3, 4 and 5 we show several spectral re- 
gions where we can see the best model synthesis 
in comparison with the synthesis using solar abun- 
dances. We also analyzed the broadened spec- 
trum of the template late type star and found 
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Fig. 3. — Best synthetic spectral fits to the UVES spectrum of the secondary star in the Cen X-4 system 
(bottom panel) and the same for a template star (properly broadened) shown for comparison (top panel). 
Synthetic spectra are computed for solar abundances (dashed line) and best fit abundance (solid line). 



the abundances of all the elements under study 
to be close to solar (see Bodaghee et al. 2003) ex- 
cept the Li abundance which was estimated to be 
loge(Li) < 0.7. 

Abundances for all the elements are listed in Ta- 
ble 2 and referred to the solar values adopted from 
Anders & Grevesse (1989). We also give errors, 
A CT , estimated from the dispersion of the elemen- 
tal abundances obtained from the best fits to the 
various spectral features. Errors associated with 
uncertainties in effective temperature, Ax off , grav- 
ity, Alogtjf, and veiling, A vc iii ng , are also listed in 
Table 2. 

The Li feature is particularly strong in our tar- 
get resulting in an abundance of loge(Li)LTE = 
3.06 ± 0.29. We have also estimated the non-LTE 
abundance correction for this element, 

Aloge(Li) = loge(Li) NLT E - loge(Li) LTE 



from the theoretical LTE and non-LTE curves of 
growth in Pavlenko & Magazzu (1996). By using 
the derived stellar parameters and LTE Li abun- 
dance, we obtain loge(Li)NLTE = 2.98 for the sec- 
ondary star in this system. 

The Lil 6708A resonance doublet may also form 
in absorption in the atmosphere of the accretion 
disk during outburst decay (Sulcimanov & Rcbolo 
1998). Such lines of Li and possibly of other ele- 
ments could arise if the disk luminosity is roughly 
0.1-0.005 the Eddington luminosity. However, 
Cen X-4 was observed at quiescence, when its X- 
ray luminosity was close to 10 -6 the Eddington 
luminosity (Asai et al. 1996). In addition, these 
lines would be extremely broadened as a conse- 
quence of the rotation profile of the disk (~ 500- 
1000 kms -1 ). Therefore, we discard that absorp- 
tion lines formed in the accretion disk have any 
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Fig. 4. — The same as in Figure 3. The spectrum of the template is not corrected for telluric lines (atm. 
H2O); therefore, these lines appear broadened in the spectrum (histogram line) and the synthetic spectrum 
does not fit (solid line). 



significant influence on the derived chemical abun- 
dances of the secondary star. 

In the following section, we will discuss these 
results in the framework of the origin and evolu- 
tionary scenario of the Cen X-4 system. 

4. Discussion 

4.1. Heavy Elements 

The abundances of heavy elements in the se- 
condary star in Cen X-4 are super solar. In Fig. 
6 these element abundances relative to iron are 
shown in comparison with the Galactic abundance 
trends of these elements in the relevant metallicity 
range, taken from Feltzing & Gustafsson (1998) 
and Bodaghcc ct al. (2003). The error bars in 
the element abundance ratios ([E/Fe], see Table 
3) take into account how individual element abun- 



dances depend on the various sources of uncer- 
tainty when dealing with abundance ratios. The 
[Ca/Fe] and [Al/Fe] ratios of the secondary are 
consistent with abundances of stars with similar 
iron content, while Ni and Ti appear to be mod- 
erately enhanced. In Table 3 we show the element 
abundance ratios in the secondary star in Cen X-4 
and average values in stars with iron content in 
the range —0.13 < [Fe/H] < 0.33, the comparison 
sample, corresponding to a la uncertainty in the 
iron abundance of the companion star. Whereas 
[Ca/Fe] and [Al/Fc] ratios are consistent with av- 
erage values, the [Ni/Fe] and [Ti/Fe] ratios are la 
higher than the average values of the stars in the 
comparison sample. 

It has been proposed that LMXBs such as 
Cen X-4 begin their lives as wide binaries and 
evolve through a common-envelope phase where 
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Fig. 5. — The same as in Figure 3. 



the companion spirals into the massive star's en- 
velope (see, for example, Van den Heuvel 1983; 
de Kool et al. 1987; Nelemans & Van den Heuvel 
2001). The helium core of the massive star contin- 
ues its evolution until dying as a supernova leaving 
a neutron star or a black hole remnant. Part of 
the ejected mass in the SN explosion could be cap- 
tured by the companion star, as has been found in 
LMXBs like Nova Sco 94 (Israelian et al. 1999) and 
A0620-00 (Gonzalez Hernandez et al. 2004). We 
have considered this possibility here in the case of 
the Cen X-4 system. 

4-1.1. Spherical SN explosion 

A binary system like Cen X-4 (with neutron 
star mass Mns = 0.5 — 2.5 M© and the com- 
panion star mass in the range M 2 = 0.04-0.58 
Mq) will survive a spherical SN explosion if the 
ejected mass, AM = M Hc -M NS < (M He +M 2 )/2. 
Therefore, initial masses for a neutron star and 



the secondary star of ~ 1.4 M© and ~ 0.8-1 M© 
respectively would imply a mass of the He core 
before the SN explosion M He < 3.6-3.8 Mq. Us- 
ing the expressions given by Portegies Zwart et 
al. (1997, and references therein), we inferred a 
mass of the progenitor star, Mi ~ 16 M©, and 
a radius Ru e ~ 0.8 Rq for the helium core. As 
in Gonzalez Hernandez et al. (2004), we can also 
estimate the amount of the ejected material in a 
spherical explosion that could be captured by the 
companion as coming from a central point since 
the pre-SN orbital separation between the binary 
components would be a <~ 1.6-2 Rq (assuming 
M He ~ 4 Mq, M ns ~ 1.5 Mq and M 2 ~ 0.5-0.8 
Mq). The explosion of such a massive primary 
would have taken place only ~ 7 x 10 6 yr after 
the formation of the system (Brunish & Truran 
1982). At that time, the radius of an 0.8 M© 
secondary star would be ~ 1.2 Rq (DAntona & 
Mazzitelli 1994). If we consider a spherically sym- 
metric supernova explosion, taking into account 
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Fig. 6. — Abundances of the secondary star in Cen X-4 (wide crosses) in comparison with the abundances 
in G and K metal-rich dwarf stars. Galactic trends of Ca, Ni, and Ti were taken from Bodaghee et al. 
(2003) while Al from Feltzing & Gustafsson (1998). The size of the cross indicates the size of error. The 
dashed-dotted lines indicate solar abundance values. 



the fraction of solid angle subtended by the com- 
panion and assuming different capture efficiency 
factors, /capture, the amount of mass deposited on 
the secondary can be computed as 

m add = AM(7ri?2/47rao)/ capturo 

In Table 4 we present the expected abundances 
in the atmosphere of the secondary star after con- 
tamination from the progenitor of the compact ob- 
ject for two plausible secondary masses (0.8 M Q 
and 0.5 M Q ) at the time of the explosion of the pri- 
mary. We have used a spherically symmetric core- 
collapse explosion model of He core of mass Mn e ~ 
4 Mq and explosion energy of e = 10 51 ergs. This 
energy is deposited promptly in the central region 
of the progenitor core to generate a strong shock 
wave. The subsequent propagation of the shock 
wave is followed through a two-dimensional hydro- 
dynamical code, as has been done for aspherical 
explosion models by Maeda et al. (2002). Nucle- 
osynthesis is solved by detailed nuclear reaction 
network calculations including 222 isotopes, from 



1 H up to 71 Gc. In the simulations, the directional 
distribution of the input energy at the mass cut 
is given in a parametric way to allow asymmetric 
explosions. One half of the energy is deposited as 
kinetic energy, and the corresponding initial veloc- 
ities (v r , v z ) are imposed as v z = az and v r = (ir 
(here r is the equatorial direction, z the jet direc- 
tion, and v r and v z the initially imposed veloc- 
ity components along each direction). Therefore, 
the initially imposed velocity is a factor of a/f3 
larger in the jet direction than in the equatorial 
direction. In the spherical symmetric model, the 
ratio a/fl is set to be unity, therefore the model 
is similar to canonical ID explosion calculations 
(see, e.g., Umeda et al. 2000). The details of our 
explosion calculations can be found in Maeda et 
al. (2002). Our present models are basically the 
same as their models, except for the progenitor's 
mass and the explosion energy. In our modeling, 
we have chosen different capture efficiency factors, 
/capture, between and 1. However, simulations of 
type la SNe suggest that the SN blast wave may 
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Table 3 

Element abundance ratios in the secondary star in Cen X-4 and in the comparison sample 



Element 


[ E / Fc ] C e„X-4 


^[E/Fe],CenX-4 


[E/Fe] atarB 


CT stars 


<^er, stars 


Al 


0.07 


0.13 


0.04 


0.03 


0.006 


Ca 


-0.02 


0.18 


-0.07 


0.04 


0.008 


Ti 


0.17 


0.14 


-0.04 


0.06 


0.010 


Ni 


0.12 


0.08 


-0.02 


0.05 


0.010 



*Errors in the element abundance ratios ([E/Fc]) in the secondary star in Cen X-4. 



Note. — [E/Fe] stars indicates the average values calculated for stars with iron con- 
tent in the range —0.13 to 0.33 corresponding to la in the [Fc/H] abundance of the 
secondary star in Cen X-4. Ca, Ti and Ni for the comparison sample have been taken 
from 29 stars in Bodaghee et al. (2003) while Al from 28 stars in Feltzing & Gustafs- 
son (1998). The uncertainty in the average value of element abundance ratios in the 
comparison sample is obtained as A CTjSt ars — Tatars/ VlV where <r a tars is the standard 
deviation of the measurements and N, the number of stars. 



induce mass-loss in the secondary instead of mat- 
ter accretion (Marietta et al. 2000) and therefore, 
high values for the capture efficiency may not be 
appropriated. We also assume that the matter 
captured is well-mixed in the convective envelope 
of the companion star. The mass of the convective 
zone of the secondary star was fixed at m conv = 0.5 
and 0.6 M Q from the evolutionary tracks for a 0.5 
and 0.8 M Q secondary star, respectively, with an 
age of - 7 x 10 6 yr (D'Antona & Mazzitelli 1994). 

The mass cut, M cut , that divides the ejecta 
from the collapsing core, is low enough for a sig- 
nificant fraction of iron to be ejected and conse- 
quently captured by the secondary star. The sec- 
ondary star could therefore have had an initial iron 
abundance lower than its current observed abun- 
dance. In order to check this possibility we have 
assumed different initial mctallicities for the sec- 
ondary and primary star in the system, including 
low mctallicities (given the high galactic latitude, 
b <~ 24° , of the system that might suggest a halo 
origin). We have taken the initial element abun- 
dances from the galactic trends (Allcndc Prieto 
et al. 2004; Stephens & Boesgaard 2002) for four 
given iron abundances. We have used nucleosyn- 
thetic products of the solar mctallicity explosion 
model since the abundance pattern of a-nuclei can 
be similar to the solar abundance (Umeda et al. 
2000). However, the odd Z to even Z element 
ratios become smaller for lower metallicity, and 
hence we have excluded 27 A1 from low mctallicity 
models. 

The expected abundances for the solar abun- 



dance model presented in Table 4 can fit very well 
with the observed abundances in the secondary 
star for lower mass cuts. However, these element 
abundances are strongly sensitive to the mass cut, 
and therefore for a slighlty higher value of this 
parameter the expected abundances cannot repro- 
duce the observed abundances. This strong depen- 
dence on the mass cut makes the low metallicity 
models unable to produce acceptable fits to the 
observed abundances. In addition, note that for 
lower metallicity models the secondary star needs 
to capture 10-20 per cent of the mass in its convec- 
tive zone from the ejecta (i.e., an /capture ~ 0.2- 
0.7) to achieve the current observed abundances. 
Whereas for the solar abundance model by cap- 
turing roughly 5% of its convective envelope mass 
the secondary star can increase sufficiently its ele- 
ment abundances (i.e., an /capture ~ 0.1-0.3). Ac- 
cording to the chemical abundance pattern in the 
secondary star, wc find as very unlikely a halo ori- 
gin. Although, strictly speaking, we may only rule 
out an origin in the metal-weak halo. The halo is 
also populated by stars that have been accreted 
from Milky Way satellites and some of these are 
of high metallicity. The only known satellites that 
are accreting onto the Milky Way and comprise a 
population as metal-rich as the Cen X-4 system 
are Sagittarius (Ibata et al. 1994, Bonifacio et al. 
2004) and Canis Major (Martin et al. 2004, Sbor- 
done et al. 2004). Sagittarius is in a polar orbit 
(Ibata ct al. 1997) and Canis Major in an equato- 
rial orbit (Martin et al. 2004). The present orbit 
of Cen X-4 (see Gonzalez Hernandez et al. 2005) is 
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Table 4 

Spherical Supernova Explosion in Cen X-4 



[E/H] EXPECTED d 











M 2 = 


-0.5 Mq 


M 2 = 


0.8 M 












.t(M ) 


M cut (M ) 


element 


[E/H] OBS. a 


[E/H] SAMPLE b 


[E/HL ° 


1.49 


1.55 


1.49 


1.55 






[Fe/H] fii 


= 










Al 


0.30 


0.25 





0.27 


0.47 


0.27 


0.48 


Ca 


0.21 


0.15 





0.23 


0.37 


0.23 


0.38 


Ti 


0.40 


0.18 





0.48 


0.19 


0.49 


0.20 


Fe 


0.23 


0.23 





0.23 


0.23 


0.23 


0.23 


Ni 


0.35 


0.20 





0.32 


0.18 


0.32 


0.19 






[Fe/H] tji = 


-0.4 










Al 


0.30 


0.25 












Ca 


0.21 


0.15 


-0.24 


0.26 


0.49 


0.27 


0.49 


Ti 


0.40 


0.18 


-0.11 


0.65 


0.26 


0.66 


0.26 


Fe 


0.23 


0.23 


-0.38 


0.23 


0.23 


0.23 


0.23 


Ni 


0.35 


0.20 


-0.31 


0.39 


0.17 


0.39 


0.17 






[Fe/H] tjj = 


-0.8 










Al 


0.30 


0.25 












Ca 


0.21 


0.15 


-0.71 


0.24 


0.50 


0.23 


0.51 


Ti 


0.40 


0.18 


-0.63 


0.67 


0.17 


0.67 


0.17 


Fe 


0.23 


0.23 


-0.78 


0.23 


0.23 


0.23 


0.23 


Ni 


0.35 


0.20 


-0.65 


0.41 


0.14 


0.41 


0.15 






[Fe/H] tjj = 


-1.4 










Al 


0.30 


0.25 












Ca 


0.21 


0.15 


-1.17 


0.23 


0.52 


0.52 


0.52 


Ti 


0.40 


0.18 


-1.16 


0.69 


0.15 


0.15 


0.16 


Fe 


0.23 


0.23 


-1.41 


0.23 


0.23 


0.23 


0.23 


Ni 


0.35 


0.20 


-1.12 


0.42 


0.13 


0.13 


0.13 



a Observed abundances of the secondary star in Cen X-4. 

b Average abundances in stars of the comparison sample (sec also Table 4). 

c Initial abundances assumed for the secondary star in Cen X-4 from Allende Prieto et al. 
(2004) and Stephens & Boesgaard (2002) 

d Expected abundances of the secondary star 



Note. — Expected abundances in the secondary atmosphere contaminated with nuclcsynthctic 
products of a 16 Mq spherically symmetric core-collapse explosion model (Mh c ~ 4 Mq) for 
two different secondary masses, M 2 . M cu t is the mass cut assumed for each model. [Fe/HL i is 
the initial mctallicity of the secondary star. The capture efficiency, /capture , has been modified in 
each model until matching the expected with the observed iron abundance. Therefore, the mass 
captured by the companion, m a dd, lies between 0.02 (/capture — 0.1) and 0.11 (/capture — 0.5) 
for the extreme values of M 2 = 0.5 Mq, M cu t = 1-49 Mq and [Fe/H] t i = 0, and M 2 = 0.8 Mq, 
Mcut = 1-55 Mq and [Fe/H] + )j = -1.4 
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not compatible with that of either of these galax- 
ies, thus it seems likely that Cen X-4 originated in 
neither of them, but a detailed kinematic study is 
needed to reject this possibility. 

4-1.2. Aspherical SN explosion 

So far we have been considering a spherically 
symmetric SN explosion, but the kinematic prop- 
erties of this system might suggest an aspherical 
SN explosion, as has been proposed for Nova Sco 
94 (Brown ct al. 2000). The maximum ejected 
mass in a spherical SN explosion, AM <~ 2.4 
M©, leads to a system velocity of ~139 km s _1 . 
Gonzalez Hernandez et al. (2005) measured the 
proper motion of the system and derived a space 
velocity of w sys > 300 kms -1 , i.e., considerably 
larger than what could be acquired in the case of 
a spherical SN explosion. In addition, the Galac- 
tic space velocity components (U, V, W) of the sys- 
tem are significantly different from the mean val- 
ues that characterize the kinematics of stars be- 
longing to the thin and the thick disks of the 
Galaxy. We conclude that the nascent neutron 
star must have received an additional kick, most 
probably at birth. Theorists have been trying to 
explain these kicks, whose physical origin remains 
unclear (see Lai et al. 2001 for a review, and ref- 
erences therein). A natal kick could be generated 
by an asymmetric explosion due to global hydro- 
dynamic perturbations in the supernova core, or 
it could be a result of asymmetric neutrino emis- 
sion in the presence of superstrong magnetic fields 
(B > 10 15 G) in the proto-neutron star and/or 
non-standard v physics. In the latter cases, an 
asymmetry of ~1% in the neutrino emission is 
sufficient to account for a kick of 300 km s _1 ; 
therefore, a non-spherically symmetric SN explo- 
sion is not necessary. Recently, some authors have 
started to explore the asymmetric-collapse kick 
mechanisms from global hydrodynamic pertuba- 
tions. Two and three dimensional studies pre- 
dict strong neutron star kicks from this mecha- 
nism (Burrows & Hayes 1996; Janka et al. 2004) 
even in excess of 1000 km/s. On the contrary, 
Fryer (2004), on the basis of three dimensional su- 
pernova simulations, found that even the most ex- 
treme asymmetric collapses do not produced final 
neutron star velocities above 200 kms -1 . There- 
fore, an aspherical SN explosion cannot be ruled 
out and might be the explanation for the high 



space velocity of the system. 

Nucleosynthesis in aspherical explosions has 
been examined by Maeda et al. (2002), who 
showed that the chemical composition of the ejecta 
is strongly dependent on direction. In partic- 
ular, the main effect of the asphericity is that 
elements produced by the strong a-rich freeze- 
out are greatly enchanced relative to iron (e.g., 
[Ti/Fe]). This effect is more pronounced in the 
jet direction, where the kinetic energy is higher. 
Therefore, assuming that the jet is perpendicular 
to the orbital plane, i.e., that the secondary star is 
located in the equatorial plane of the helium star, 
elements like Ti, Ni, and Fe are mainly ejected in 
the jet direction, while Al, O, Si, S, and Mg are 
preferentially ejected near the equatorial plane. 
Consequently, we should have found a high Al 
abundance in the companion star in comparison 
with the other elements studied but this is not 
the case. To compare these effects with our ob- 
servations of Cen X-4 quantitatively, we compute 
an aspherical explosion of a 4M Q He core with 
the explosion energy 10 51 ergs. The asphericity is 
induced by the parametric way (see §4.1.1) with 
the asphericity parameter a/ [3 — 2, i.e., the initial 
shock velocity at the mass cut is 2 times larger 
along the jet direction than in the equatorial di- 
rection. In table 5 we show the expected clement 
abundances in the secondary star after pollution 
from the ejecta of an aspherical SN explosion of a 
4 Mq He core. Contrary to the spherical case, Al 
is dramatically enhanced in all the model compu- 
tations and this effect is even stronger for higher 
mass cuts. Therefore, for the solar abundance 
model, none of the model computations produces 
acceptable fits, even for lower mass cuts. Note 
that explosive burning in aspherical SN is weak, 
ejecting little Ti and Ni, in the equatorial direc- 
tion. Thus, the mass cut appropiate to the system, 
i.e., ~ 1.4 Mq, is not deep enough to eject these 
elements to explain the large values of [(Ti, Ni) /H] 
observed. The mass cut in the aspherical case is 
placed at relatively weak burning region as com- 
pared with the spherical model. Therefore, the 
dependence of the model abundance on changing 
the mass cut is different from the spherical model, 
which also makes the aspherical model difficult to 
fit the observed abundances. 

Note that all the above results of the aspherical 
explosion rely on the assumption that the jet di- 
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Table 5 

Aspherical Supernova Explosion in Cen X-4 



[E/H] EXPECTED d 











M2 — 


-0.5 Mg 


M 2 = 


0.8 M 












it(M ) 


M cut (M Q ) 


ELEMENT 


[E/H] OBS. a 


[E/H] SAMPLE 15 


[E/H] j ; c 


1.49 


1.55 


1.49 


1.55 






[Fe/H] tij 


= 










Al 


0.30 


0.25 





0.59 


0.84 


0.59 


0.84 


Ca 


0.21 


0.15 





0.43 


0.24 


0.43 


0.24 


Ti 


0.40 


0.18 





0.21 


0.20 


0.21 


0.20 


Fo 


0.23 


0.23 





0.23 


0.23 


0.23 


0.23 


Ni 


0.35 


0.20 





0.18 


0.22 


0.18 


0.22 






[Fe/H] tji = 


-0.4 










Al 


0.30 


0.25 












Ca 


0.21 


0.15 


-0.24 


0.56 


0.28 


0.56 


0.28 


Ti 


0.40 


0.18 


-0.11 


0.28 


0.26 


0.28 


0.26 


Fe 


0.23 


0.23 


-0.38 


0.23 


0.23 


0.23 


0.23 


Ni 


0.35 


0.20 


-0.31 


0.15 


0.22 


0.15 


0.22 






[Fe/H] tjj = 


-0.8 










Al 


0.30 


0.25 












Ca 


0.21 


0.15 


-0.71 


0.59 


0.25 


0.59 


0.25 


Ti 


0.40 


0.18 


-0.63 


0.20 


0.17 


0.20 


0.17 


Fc 


0.23 


0.23 


-0.78 


0.23 


0.23 


0.23 


0.23 


Ni 


0.35 


0.20 


-0.65 


0.12 


0.21 


0.12 


0.21 






[Fe/H] tjj = 


-1.4 










Al 


0.30 


0.25 












Ca 


0.21 


0.15 


-1.17 


0.61 


0.25 


0.61 


0.25 


Ti 


0.40 


0.18 


-1.16 


0.19 


0.15 


0.19 


0.15 


Fe 


0.23 


0.23 


-1.41 


0.23 


0.23 


0.23 


0.23 


Ni 


0.35 


0.20 


-1.12 


0.11 


0.20 


0.11 


0.20 



a Observed abundances of the secondary star in Cen X-4. 

b Average abundances in stars of the comparison sample (sec also Table 4). 

c Initial abundances assumed for the secondary star in Cen X-4 from Allende Prieto et al. 
(2004) and Stephens & Boesgaard (2002) 

d Expected abundances of the secondary star 



Note. — Expected abundances in the secondary atmosphere contaminated with nuclcsynthctic 
products of a 16 Mq non- spherically symmetric core-collapse explosion model (Mh c ~ 4 Mq) for 
two different secondary masses, M%. M cut is the mass cut assumed for each model. [Fe/H]^. i is 
the initial mctallicity of the secondary star. The capture efficiency, /capture , has been modified in 
each model until matching that expected with the observed iron abundance. Therefore, the mass 
captured by the companion, m a dd, lies between 0.04 (/capture — 0.2) and 0.22 (/capture — 0.8) 
for the extreme values of M 2 = 0.5 Mq , M cu t = 1-49 Mq and [Fe/H] + i = 0, and M 2 = 0.8 Mq, 
Mcut = 1.55 Mq and [Fe/H] + >j - -1.4 
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rection is perpendicular to the orbital plane, there- 
fore to the equatorial direction of the progenitor 
star since the star may well co-rotate with the bi- 
nary orbit. The assumption that the jet is almost 
perpendicular to the equatorial direction is justi- 
fied because recent core-collapse theories suggest 
the jets, if formed, propagate along the progen- 
itor's rotational axis (e.g., Fryer & Warren 2004 
and references therein). This is also supported 
by a Hubble Space Telescope image of SN 1987A 
(Wang et al. 2002). The ejecta of SN 1987A is 
elongated toward the direction perpendicular to 
the innermost circumstellar ring, which is proba- 
bly along the equatorial direction (e.g., Chevalier 
& Dwarkadas 1995). 

In the case of low initial metallicity of the 
secondary star, higher values for the capture ef- 
ficiency (i.e., /capture ~ 0.3-1) are required to 
achieve the observed Fe abundance since this ele- 
ment is mainly ejected perpendicularly to the or- 
bital plane. Despite this requirement, these sim- 
ulations cannot fit the observed abundances even 
for high mass cuts. 

In summary, a spherically symmetric SN explo- 
sion of a 4 Mq He core leaving an NS gives a better 
fit than an asphcrical model. A relatively low pri- 
mary star (Mi < 20 Mq) seems to leave an NS 
and to explode as a normal SN (with an energy 
of e = 10 51 ergs, e.g., Nomoto et al. 2003). Our 
result would suggest that a normal SN would be 
less aspherical than a hypernova, as proposed for 
the black hole in the LMXB Nova Sco 94 (Brown 
et al. 2000). However, we cannot rule out that 
the companion star did not capture a significant 
amount of matter from the ejecta of a SN explo- 
sion, and therefore that the observed abundances 
correspond to its primordial abundances. 

4.2. Li abundance 

The companion star has an unexpectedly high 
7 Li abundance close to the cosmic value in the 
Galactic disk (Martin et al. 1994a), substantially 
higher than field late-type main sequence stars of 
similar mass. Convective mixing during pre-main 
sequence and main sequence evolution is expected 
to produce significant lithium depletion in the at- 
mospheres of such stars (see Fig. 7), thus a possi- 
ble explanation for the Li over-abundance in the 
secondary is that the Cen X-4 system is younger 
than the Pleiades cluster whose age is ~ 120 Myr 
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Fig. 7. — Li abundance of the secondary star 
in Cen X-4 (filled circle) in comparison with the 
abundances of different rotating Pleiades dwarf 
stars versus effective temperature from Garcia 
Lopez et al. (1994). The sizes of the circles are re- 
lated to wsini. Stellar masses have been assigned 
following the evolutionary tracks of D'Antona & 
Mazzitelli (1994) for 10 8 yr. 



(Martin et al. 1998; Stauffer et al. 1998). One 
could estimate an upper limit to the age of the sys- 
tem under the assumption that there is no mech- 
anism able to enrich the atmosphere of the sec- 
ondary star with freshly synthesized Li nuclei. In 
Fig. 8 we show the age of the system for different 
mass transfer rates and initial secondary masses by 
following the Li depletion in low mass stars accord- 
ing to evolutionary tracks of Siess et al. (2000). 
We have adopted the NLTE Li abundance and the 
present value for the initial metallicity of the sys- 
tem since a possible metallicity enrichment due 
to a SN explosion would have taken place early 
enough (at ~7 Myr) and it would make the sec- 
ondary follow the Li depletion as a metal-rich star. 
In any case, this metallicity effect is only rele- 
vant for masses higher than 0.8 Mq. Such high 
Li abundance restricts the age of the system to 
the range ~ 10-70 Myr. This age is quite short 
compared to the age assumed in the proposed evo- 
lutionary scenarios for the system, where an age 
of a few times 10 9 yr is required (e.g. Chevalier 
et al. 1989). These age constrains are obtained 
using models that do not consider the impact of 
high rotation on the Li depletion. According to 
Martin & Claret (1996) high rotation may inhibit 
the rate of Li depletion; therefore, the constraints 
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on age should be strictly considered lower 
limit. The effect of rotation on Li depletion has 
been measured in Pleiades stars with masses in the 
range 0.6-0.8 M (see, for example, Garcia Lopez 
et al. 1994). The kinematic properties of Cen X-4 
suggest that it could have been at the inner re- 
gions of the Galactic disk about 100-200 Myr ago 
(Gonzalez Hernandez et al. 2005). If the neutron 
star formed at that time as a result of an SN ex- 
plosion, this might explain the high Li abundance 
of the system without postulating a Li production 
mechanism. 

For such a young age the effect of the mass 
transfer in the secondary's mass is negligible, and 
initial secondary masses higher than 0.8 M can 
therefore be ruled out. Note that the mass trans- 
fer rate inferred from the observed quiescent X-ray 
luminosity (~ 2.4 x 10 32 ergs" 1 , Asai et al. 1996) 
is very low ~ 10~ 6 the Eddington rate (~ 3 x 10~ 9 
M Q yr" 1 , Yi & Narayan 1997). On the other 
hand, Menou & McClintock (2001) have suggested 
a mass transfer rate of ^3 x 10~ 10 M Q yr -1 based 
on optical-UV emission spectrum of the system. 
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Fig. 8. — Upper limit to the age of the Cen X-4 
system according to the Li abundance of the sec- 
ondary star versus initial stellar masses of the 
secondary star, assuming different mass transfer 
rates. Li depletion has been followed according to 
the evolutionary tracks of Siess et al. (2000). 

High Li abundances appear to be a common 
property of late-type companions in Soft X-ray 
Transients (SXTs), characterized by strong X-ray 
outbursts (ussually lasting for several months) fol- 
lowed by long quiescent periods of ^10-50 yr. 
Different models have been proposed to explain 



these Li over- abundances in terms of spallation 
reactions or a-a reactions. Martin et al. (1994b) 
proposed that Li production could come from nu- 
clear reactions in the accretion disk during out- 
burst and being deposited on the companion via a 
disk wind. They argued that the observed abun- 
dance must correspond closely to the long term 
average since the depletion time scale {td > 10 7 
yr) is much longer than the recurrence time of out- 
burst (~ 30 yr). Later on, Yi & Narayan (1997) 
investigated the same scenario, in an advection- 
dominated accretion flow (ADAF), rather than a 
Shakura-Sunyacv disk not only during outburst 
but also in quiescence, although Li production 
during outburst may even dominate over the pro- 
duction in quiescence. This model can explain 
the observed Li abundance in Cen X-4 only if the 
mass-accretion rate in quiescence is ~ 10~ 3 times 
the Eddington rate and if there is a "propeller" ef- 
fect. More recently, Guessoum & Kazanas (1999) 
have proposed Li production via neutron spal- 
lation of CNO elements in the secondary star. 
In this scenario, a flux of neutrons produced in 
the hot accretion flow around the compact object 
could escape easily and impinge on the secondary 
surface and then produce spallation of heavier el- 
ements. A short exposure to the above flux, i.e., 
during outbursts, would suffice to produce the ob- 
served enhancement. 

A main sequence K star of mass ~0. 1-0.8 M Q 
would destroy all its primordial Li abundance in 
< 80 Myr. Yi & Narayan (1997, see also references 
therein) have discussed in some detail the question 
of Li depletion in the atmosphere of the companion 
star once it has been deposited/produced there, 
considering the possibility of a steady state, in 
which the enrichment and depletion rates should 
be equal, as suggested by Martin et al. (1994b). 
Martin et al. (1994a) remarked that isotopic ra- 
tios 7 Li/ 6 Li as low as 5 are expected and hence 
would provide a proof of the spallation scenario. 
Unfortunately, at such high Li abundances there 
is degeneracy between Li abundance and isotopic 
ratio and consequently, different Li abundances 
and different isotopic ratios give rise to the same 
line strength. Higher S/N spectra are required to 
extract this information from the 6708 A Li line 
although the high rotational velocity of the star 
might make it very difficult. In addition, Gues- 
soum & Kazanas (1999) suggested the possible en- 
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hancement of other light elements such as Be and 
B that might be analyzed using UV transitions. 

5. Conclusions 

We have obtained a high quality spectrum of 
the secondary star in Cen X-4 and derived atmo- 
spheric parameters and chemical abundances of 
several elements. We have applied a technique 
that provides a determination of the stellar pa- 
rameters taking into account any possible veiling 
from the accretion disk. We find T e s = 4500 ± 100 
K, log g = 3.9 ± 0.3, and a disk veiling (defined 
as -Fdisk/-Ftotai) of less than 55 per cent at 5000 A 
and decreasing towards longer wavelengths. 

The abundances of Fe, Ca, Ti, Al, and Ni are 
super solar. The abundance ratios of each element 
with respect to Fe were compared with these ratios 
in late-type main sequence metal-rich stars. Mod- 
erate anomalies for Ni and Ti have been found. A 
comparison with element yields from spherically 
symmetric supernova explosion models suggests 
that the secondary star could have captured part 
of the ejecta from a supernova that also originated 
the compact object in Cen X-4 if the primordial 
abundances of the secondary star were solar. The 
observed abundances can be explained if a progen- 
itor with a ~ 4 Mq helium core exploded with a 
mass cut ~ 1.5 Mq, such that a significant amount 
of iron could escape from the collapse of the inner 
layers, and hence the ejecta could also have en- 
hanced the iron abundance in the secondary star. 

The high galactic latitude of the system might 
suggest that the system could have originated 
in the Galactic halo and later enriched by nu- 
cleosynthetic products ejected in the supernova 
but our model computations for low metallicities 
show that this possibility is unlikely. In addition, 
the kinematic properties of this system (Gonzalez 
Hernandez et al. 2005) suggest that a natal kick 
was imparted to the neutron star at birth due to 
an asymmetric mass ejection and/or an asymme- 
try in the neutrino emission. We have also in- 
spected a non-spherically symmetric SN explosion 
model but this model cannot produce acceptable 
fits to the observed abundances in the secondary 
star. 

The Li abundance in the secondary star in 
Cen X-4 is dramatically high in comparison with 
field late-type main sequence stars, possibly indi- 



cating that this is a young system (~10-80 Myr). 
Alternatively, if the system is much older a Li pro- 
duction mechanism is needed to balance the large 
destruction of Li expected in such low-mass stars. 

6. Acknowledgements 

This work has made use of the VALD database 
and IRAF facilities and has been partially financed 
by the Spanish Ministry projects AYA2001-1657, 
AYA2002-03570 and AYA2002-0036. 

REFERENCES 

Anders, E.,& Grevesse, N. 1989, Geochim. Cos- 
mochim. Acta, 53, 197 

Asai, K., Dotani, T., Mitsuda, K., Hoshi, R., 
Vaughan, B., Tanaka, Y., & Inoue, H. 1996, 
PASJ, 48, 257 

Allende Prieto, C, Barklem, P. S., Lambert, D. 
L., & Cunha, K. 2004, A&A, 420, 183 

Bodaghee, A., Santos, N. C, Israelian, G., & 
Mayor, M. 2003, A&A, 404, 715 

Bonifacio, P., Sbordonc, L., Marconi, G., 
Pasquini, L., & Hill, V. 2004, A&A, 414, 503 

Brown, G. E., Lee, C.-H, & Bethe, H. A. 1999, 
New A, 4, 313 

Brunish, W. M.,& Truran, J. W. 1982, ApJS, 49, 
447 

Burrows, A., & Hayes, J. 1996, Phys. Rev. Lett., 
76, 352 

Canizares, C. R., McClintock, J. E., & Grindlay, 
J. E. 1980, ApJ, 236, 55 

Chevalier, C, Ilovaisky, S. A., van Paradijs, J., 
Pedersen, H., & van der Klis, M. 1989, A&A, 
210, 114 

Chevalier, R. A., & Dwarkadas, V. V. 1995, ApJ, 
452, L45 

Conner, J. P., Evans, W. D., & Belian, R. D. 1969, 
ApJ, 157, 157 

D'Antona, F.,& Mazzitelli, I. 1994, ApJS, 90, 467 

Fcltzing, S.,& Gustafsson, B. 1998, A&AS, 129, 
237 



16 



Fryer, C. L. 2004, ApJ, 601, L175 

Fryer, C. L., & Warren, M. S. 2004, ApJ, 601, 391 

Garcia Lopez, A., Rebolo, R., & Martin, E. L. 
1994, A&A, 282, 518 

Gessoum, N., & Kazanas, D. 1999, ApJ, 512, 332 

Gonzalez Hernandez, J. I., Rebolo, R., Israelian, 
G., Casares, J., Maeder, A., & Meynct, G. 2004, 
ApJ, 609, 988 

Gonzalez Hernandez, J. I., Rebolo, R., 
Penarrubia, J., Casares, J., & Israelian, 
G. 2005, submitted. 

Ibata, R. A., Wysc, R. F. G., Gilmorc, G., Irwin, 
M. J., & Suntzeff, N. B. 1997, AJ, 113, 634 

Ibata, R. A., Gilmore, G., & Irwin, M. J. 1994, 
Nature, 370, 194 

Israelian, G., Rebolo, R., Basri, G., Casares, J., & 
Martin, E. L. 1999, Nature, 401, 142 

Janka, H.-Th., Schcck, L., Kifonidis, K., Midler, 
E., & Plcwa, T. 2004, astro-ph/0408439 

Kaluzienski, L. J., Holt, S. S., & Swank, J. H. 
1980, ApJ, 241, 779 

de Kool, M., van den Heuvel, E. P. J., & Pylyser, 
E. 1987, A&A, 183, 47 

Kurucz, R. L., Furenild, I., Brault, J., & Tester- 
man, L. 1984, Solar Flux Atlas from 296 to 1300 
nm, NO AO Atlas 1 (Cambridge: Harvard Univ. 
Press) 

Lai, D., Chernoff, D. F., & Cordes, J. M. 2001, 
ApJ, 549, 1111 

Maeda, K., Nakamura, T., Nomoto, K., Mazzali, 
P. A., Patat, F., & Hachisu, I. 2002, ApJ, 565, 
405 

Marietta, E., Burrows, A., & Fryxell, B. 2000, 
ApJS, 128, 615 

Martin, N. F., Ibata, R. A., Bellazzini, M., Irwin, 
M. J., Lewis, G. F., & Dehnen, W. 2004, MN- 
RAS, 348, 12 

Martin, E. L., Rebolo, R., Casares, J., & Charles, 
P. A. 1994a, ApJ, 435, 791 



Martin, E. L., Spruit, H. C, & van Paradijs, J. 
1994b, A&A, 291, L43-L46 

Martin, & Claret, A, 1996, A&A, 306, 408 

Martin, E. L., Basri, G., Gallegos, J. E., Rebolo, 
R., Zapatero-Osorio, M. R., & Bejar, V. J. S. 
1998, ApJ, 499, 61 

Matsuoka, M., Inoue, H., Koyama, K., Mak- 
ishima, K., Murakami, T., Oda, M., Ogawara, 
Y., Ohashi, T., Shibazaki, N., Tanaka, Y., & 10 
coauthors 1980, ApJ, 240, 137 

McClintock, J. E., & Remillard, R. A. 1990, ApJ, 
350, 386 

Menou, K.,& McClintock, J. E. 2001, A&A, 557, 
304 

Nelemans, G., & Van den Heuvel, E. P. J. 2001, 
A&A, 376, 950 

Nomoto, K, Maeda, K, Umeda, H, Ohkubo, T, 
Deng, J, & Mazzali, P. 2003, A Massive Star 
Odyssey: From Main Sequence to Supernova, 
Proceedings of IAU Symposium #212, held 
24-28 June 2001 in Lanzarote, Canary island, 
Spain. Edited by Karel van der Hucht, Artemio 
Herrero, and Esteban, Csar. San Francisco: As- 
tronomical Society of the Pacific, p. 395 (astro- 
ph/0209064) 

van Paradijs, J., Vcrbunt, F., Van der Linden, T., 
Pedersen, H., & Wamsteker, W. 1980, ApJ, 241, 
161 

Pavlenko, Ya. V., & Magazzu, A. 1996, A&A, 311, 
961 

Piskunov, N. E., Kupka, F., Ryabchikova, T. A., 
Weiss, W. W., & Jeffery, C. S. 1995, A&AS, 
112, 525 

Podsiadlowski, P., Nomoto, K., Maeda, K., Naka- 
mura, T., Mazzali, P., & Schmidt, B. 2002, 
ApJ, 567, 491 

Portegies Zwart, S. M., Verbunt, F., & Ergma, E. 
1997, A&A, 321, 207 

Sbordonc, L., Bonifacio, P., Marconi, G., Zaggia, 
S., & Buonanno, R. 2004, A&ALetters, in press 
(astro-ph/0412098) 



17 



Shahbaz, T., Naylor, T., & Charles, P. A. 1993, 
MNRAS, 265, 655 

Siess, L., Dufour, E., & Forestini, M. 2000, A&A, 
358, 593 

Sncdcn, C. 1973, PhD Dissertation (Univ. of Texas 
at Austin) 

Stauffcr, J. R., Schultz, G., & Kirkpatrick, J. D. 
1998, ApJ, 499, 199 

Stephens, A., & Boesgaard, A. M. 2002, ApJ, 123, 
1647 

Suleimanov, V., & Rebolo, R. 1998, A&A, 330, 
999 

Torres, M. A. P., Casares, J., Martmez-Pais, & 
Charles, P. A. 2002, MNRAS, 334, 233 

Umeda, H., Nomoto, K., Nakamura, T. 2000, 
The First Stars. Proceedings of the MPA/ESO 
Workshop held at Garching, Germany, 4-6 Au- 
gust 1999. Achim Weiss, Tom G. Abel, Vanessa 
Hill (eds.). Springer, (astro-ph/9912248) 

Van den Hcuvel, E. P. J. 1983, in Accretion-driven 
Stellar X-ray Sources, eds. W. H. G. Lewin & 
E. P. J. Van den Heuvel, Cambridge University 
Press, p. 303 

Wang, L. et al. 2002, ApJ, 579, 671 

Yi, I. & Narayan, R. 1997, ApJ, 486, 363 



This 2-column preprint was prepared with the AAS IATjrjX 
macros v5.2. 



18 



